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C^) ' Properties of QCD at finite chemical potential (fi) are extracted using Taylor series ex- 

pansions. The continuum limits of lattice results are presented. The result of expanding the 
free energy density, i. e., the pressure, to 6th order in the expansion is shown. The Taylor 
coefficients of the chiral condensate are also shown. Relations between various Taylor coef- 
q ' ficients are demonstrated. All this information is utilised to remove various lattice artifacts 

, from the determination of the Wroblewski parameter in strangeness production. 
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Lattice simulations of QCD are possible whenever the determinant of the Dirac 
<^ ' operator, M, is positive. The proof of positivity amounts to constructing an operator 

Q such that QMQ" 1 = M'. Once such a Q has been obtained, the result that Det M 
is real follows trivially. For QCD at zero chemical potential (// = 0), Q = 75. In 
the limit when isospin symmetry is exact, and an isovector chemical potential ^3 is 
switched on, an isospin flip is the appropriate Q, thus showing that QCD at finite /Z3 
is amenable to direct lattice simulation.^ A chiral model analysis found the phase 
diagram in the T-fi^ plane P and lattice simulations^ verified these results soon 
thereafter. At finite baryon chemical potential (p) there is no such operator, and 
direct simulations are impossible. CP symmetry nevertheless dictates that the free 
energy remains real, and that normal thermodynamics is obtained^ 

Many methods have been developed to explore the interesting physics in the 
T-\i planed Here I discuss high-order Taylor expansions of the pressure and first 
results of systematic Taylor expansions of many other quantities. The first results 
on the Taylor expansion of the free energy (pressure) were presented irP and that 
of meson correlators inP 1 This report is in two section — the first deals with the 
Taylor expansion of the pressure and the second with that of the quark condensate 
and related quantities. Following these two sections is a brief summary. 

§2. The pressure 

jJEJJHl th e following series expansion is utilised for the pressure, 

P(T, (J,) = P(T, 0) + X uu(T)fl 2 + -^Xuuuu(T)fl A + ^rXuuuuuu(T)fl 6 + ■■■ (2-1) 

where the pressure, P of a system at temperature T is given in terms of the free 
energy, F, by P = —F/V, where V is the volume. A similar expansion is used by 
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the Bielefeld-Swansea group CP symmetry forces the odd terms in the expansion 
to vanishP The neglected terms are meant to be of higher order in /i. The Taylor 
coefficients are the generalised susceptibilities, which are best defined by introducing 
one chemical potential for each flavour of quarks. This is allowed since flavour is a 
conserved quantum number in QCD. Then the susceptibilities of order n are defined 
as the partial derivatives 

1 d n F 

" ~V d Hl ■ ■ ■ dfAf, 

where {fi, ■ ■ ■ f n } are flavour labels. It is then clear that two assumptions have been 
made in writing eq. (|2-1[) . first that we consider the isospin symmetric case with 2 
flavours of quarks, and second that only the diagonal susceptibilities (/i = f% = 
■ ■ ■ = f n ) need to be retained. Both these restrictions are non-essential, but they are 
justified by the respective observations that flavour asymmetry plays a negligible role 
in the problem^) and that the off-diagonal susceptibilities are negligibly small. 81 ' 91 

We note here that the susceptibilities are not merely of formal interest as expan- 
sion coefficients of the pressure, but also form observables in their own right. They 
control fluctuations and the strangeness production rate^J in heavy-ion collisions. 

The Taylor expansion can be cast into a form which makes clear how good the 
approximation is, and when it breaks down. Define [i\ to be that value of fi at which 
the i-th term is equal to the i + 2-nd (for example, [i\ = \/l2x U u/Xuuuu)- Then the 
expansion in eq. (|2-1|) can be written as 

(2-3) 

where AP(T,fi) = P(T,/j,) — P(T, 0). This is manifestly well-behaved if we have 
fj>2 > fJ>% > ■ ■ ■ , and equally well-behaved series of approximations arise by neglecting 
terms of order i and above when /i <C fi*. Thus term by term improvement of the 
series is possible. Another advantage is that each susceptibility (and hence /J,*) is 
computed at zero chemical potential, whereby the continuum limit can be (and has 
been) obtained by the usual techniques of lattice gauge theoryP^ 

The other nice point about the Taylor series expansion emerges when one ex- 
amines its failure. If the sequence fi* tends to a limit then clearly the series fails 
to converge for fi = fi*. Therefore any finite //* is an estimate of the nearest phase 
boundary. Other series extrapolation methods can also be used — for example, Pade 
analysis. Results using susceptibilities up to 8th order for QCD with 2 flavours of 
light dynamical quarks in the continuum limit will be presented elsewhere. 

Like all other lattice methods in current use for this problem, the Taylor series 
expansion is limited up to the transition point (or line) nearest to fj, = 0. I be- 
lieve that it has some advantage over other methods. The primary problem with 
direct lattice computations at finite chemical potential is to gain control over phase 
fluctuations of the determinant In the reweighting methods this is achieved by 
simultaneous movement in both T and /uP^» If one could develop a method to sum 
over CP orbits in configuration space, then the sign problem would be automatically 
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Fig. 1. Continuum limits of the diagonal susceptibilities up to 6th order at T = 2T C . The bands 
enclose 66% confidence limits, x^muu"" vanishes at the 99% confidence level. The continuum 
pressure in high temperature QCD for various values of fi/T c marked (SPS corresponds to 
lx/T c = 0.45 and RHIC to 0.15). 



solvedP In the Taylor series expansion, one starts from the free energy, where the 
sum over CP orbits has already been performed to yield a real quantity. Thus, an 
expansion in the single variable fi completely captures the physics. The mechanics 
of the proof involves taking the double Taylor series expansion of P and checking 
that it gives no advantage. 

In Figure ^ we display the continuum limits of the diagonal susceptibilities and 
the EOS. The latter is compared with the Taylor series truncated at the second order, 
and with results obtained through reweighting on lattices with Nt = 4 which were 
extrapolated to the continuunP^ using the known ratio Xuu(cont) / 'xuu{Nt = 4)P^ 
These comparisons show the interesting fact that in the high temperature phase of 
the QCD plasma the computation of the linear susceptibility Xuu is sufficient to 
determine the continuum value of the pressure at fi of interest to heavy-ion experi- 
ments. This is also evidence that extrapolation in one variable of P is equivalent to 
reweighting in two variables. 

Not all observable physics in the continuum limit can be extracted by a simple 
rescaling by the factor x uu (cont) / Xuui^t = 4). At finite chemical potential lattice 
artifacts are large and connected to the fact that there is an infinity of equivalent 
prescriptions for putting chemical potential on the latticed This results in an 
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ambiguity in the definition of susceptibilities, and in the critical end point, which 
vanishes only in the continuum limit. We have estimated that this ambiguity may 
move estimates of the critical end point by an amount much larger than the present 
statistical errors.^ Hence, taking the continuum limit becomes crucial to obtaining 
this unique signature of QCD dynamics. 

The diagonal linear susceptibilities are in good agreement with perturbative com- 
putations E3>ED Consequently, AP(T, /i) is also in good agreement with perturbative 
results. 21 ' This is interesting since P(T, 0) is badly reproduced by perturbation the- 
ory. Perturbative computations of the fourth order susceptibility are not as good, 
and the off-diagonal susceptibilities fare poorly. It has been suggested that the large 
Nf behaviour of perturbation theory can be tested®''^ to throw more light on this. 



§3. The quark condensate and meson correlators 



The success of the Taylor expansion of the pressure leads us to believe that other 
quantities of interest can also be successfully investigated at finite chemical potential 
through Taylor expansions. We use the quark condensate as an example. This is 
defined to be 

C(T, n) = {W)t* = z = ~ W (3-1) 

where Z is the partition function, and we have introduced a non-standard notation, 
C, for the quark condensate which is usually denoted by (tpip). We develop this in a 
Taylor series expansion around /x = 0, 

C(T, n) = C(T, 0) + ci (T)/i + l -c 2 (T) M 2 + • • • (3-2) 

Interesting results follow immediately. 

The first derivative at a general value of fi is 

Cl(T , rt = ^l^, (3 .3) 
o[i 1 am 

where n(T, fi) is the quark number density, obtained by taking the first derivative 
of the pressure with respect to [i. We have related two seemingly different physical 
objects by interchanging orders of derivatives; this is a prototype of a Maxwell rela- 
tion in thermodynamics. Further, the Taylor expansion of the pressure shows that 
n(T,fj,) « 2xuu(T)/j,, and hence 

frr \ 2 ^ dxuu(T) 
C1 (T^) = -^-— + ••• (3-4) 

where the neglected terms start at order /i 3 . From this it immediately follows that 
ci in eq. ()3-2|) vanishes. The outline of a proof that the Taylor expansion of C in eq. 
(|3-2j) is even is clear from this. Furthermore, it also follows that the limiting value 
of ci/fi is proportional to the slope of Xuu with mass. 
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A Maxwell relation for the second Taylor coefficient C2 follows in a very similar 
manner, 

/rp\ 1 9Xuu(T) , v 

C2{T) = - — . (3-5) 

1 am 

This is of great use in lattice computations of strangeness production. It has been ar- 
gued that the Wroblewski parameter which can be extracted from heavy-ion collision 
experiments on strangeness production can be written afl^ A S (T) = Xss{T)/xuu{T). 
The strange and light quark masses are not direct physical observables, but are free 
parameters of QCD which are obtained by fitting the spectrum of mesons at T = 0. 
At present it is hard to perform a lattice computation at realistic values of the light 
quark masses. One can then compute X s at a value of the light quark mass for which 
the computation is feasible and use a Taylor expansion in the quark mass to extrap- 
olate to the physical quark mass value. The leading term in this Taylor expansion 
is then -Tc 2 /Xuu- 




Fig. 2. The renormalised second derivative of the quark condensate, c£ for T = 1.5Tc (circles), 2T C 
(boxes) and 3T C (pentagons). The continuum limit (1/-/V t 2 = 0) is consistent with zero at the 
99% confidence level at all T. 

One subtlety in these arguments is the effect of renormalisation. Since /x appears 
as a coupling to a conserved charge, it is not renormalised, and hence renormalisa- 
tion affects the computations of the susceptibilities only in so far as they need to be 
extrapolated to the continuum limit. The situation is different for the quark conden- 
sate and its Taylor coefficients. The mass is renormalised, by a multiplicative factor 
for staggered quarks. This requires a multiplicative renormalisation of every term in 
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the series in eq. (|3-2[) . This is straightforward. Since the renormalisation constant 
depends only on the ultraviolet cutoff, i. e., the lattice spacing a, and not on T and /x, 
renormalisation can be accounted for simply by making a Taylor series expansion of 
C(T, (i)/C(0, 0). The series remains even, and results for C2(T) are shown in Figure 
121 The renormalised value of c 2 is consistent with zero at the 99% confidence level. 

The Taylor expansion can also be performed in the isovector chemical potential, 
/X3. A comparison of the corresponding Taylor coefficients reveals that numerically C2 
is very nearly equal to the corresponding coefficient in the series in /X3 (in contrast, 
exact equality is obtained in a random matrix modeP^). This turns out to be a 
fairly general feature, being reproduced also in the Taylor expansions of the meson 
correlators. It is due to the observed fact that off-diagonal susceptibilities like Xud 
are small. 

We turn now to correlators of quark bilinears, which are called meson correlators 
in extension of their zero-temperature meaning. Taylor expansions of the correlation 
lengths have been computed earlier using an expression in which each correlator is 
saturated by a single mass. 7 ' We make instead a Taylor expansion of the meson 
susceptibilities. In the thermodynamic limit of the low-temperature phase of QCD 
where mesons are physical degrees of freedom, these susceptibilities are proportional 
to the inverse square of the mass.^ 

An interesting point is the fact that these meson susceptibilities can be consid- 
ered either as a sum over the temporal correlator at all temporal separation, or as 
a screening correlatoJ^ summed over all spatial separation. Above T c the effective 
representations of these two correlators are different.^ As a result, the meson sus- 
ceptibilities can also be used to extract the dependence on /j, of both masses and 
screening lengths without any need to use a spectral function in intermediate steps. 

Further relations between coefficients of different Taylor series are given by chiral 
Ward identities such as C(T, /i) = mxps(T, n), where Xps(T, /j.) is the pseudo-scalar 
susceptibility. This particular identity is quite interesting, as can be seen by writing 
out the Taylor expansion of xps — 

Xps(T, n) = Xps(T, 0) + Xps(7> + \xps(T)v 2 + • • • • (3-6) 

The chiral Ward identity then allows us to use the arguments already given for the 
series expansion of C(T, //) to argue that the series for xps is even. 
Furthermore, we obtain the interesting relations 

XPS (T) = I C2 (T) = -±- (3-7) 
m ml am 

Recall that for T > T c , the infrared cutoff is the Matsubara frequency, fl = ttT, 
whenever m < Q. As a result, the meson screening masses are given by 2fi (or its 
lattice equivalent) j2DJ23 anc i i s independent of m. The susceptibility Xuu is equal to 
the vector meson susceptibility (for staggered quarks, the identity is with a one-link 
separated meson). By the previous argument, this must be independent of m for m < 
17. This is also seen through explicit computation. 29 ' This explains the observation 
that both the quark condensate and the PS susceptibility are independent of /Up^ 
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at least to quadratic order in /i. By the connection established earlier between 
the variation of X s and C2, we therefore conclude that the Wroblewski parameter is 
insensitive to the light quark mass, provided that m < Q. At the same time, if the 
dependence of A s on the strange quark mass, m s , is considered, then, since m s m T c , 
there should be strong dependence of X s on m s . 



§4. Summary 



We have explored the physics of QCD at finite chemical potential through sys- 
tematic Taylor expansions of various quantities. This allows us to use standard 
lattice techniques to tackle this problem, and obtain the continuum limit of every 
observable. The Taylor coefficients themselves contain interesting physics, and we 
have given some examples. There are numerous relations between them — we have 
used a chiral Ward identity and a Maxwell relation as illustrations. In this report we 
have shown how the quark mass dependence of the Wroblewski parameter is related 
to physics at finite chemical potential through such relations. 
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